Физико-математические науки 








УДК 539.3: 534.11 
Т.А. АНИКИНА, И.В. БОГАЧЕВ, А.О. ВАТУЛЬЯН 
ОБ ОПРЕДЕЛЕНИИ НЕОДНОРОДНЫХ РЕОЛОГИЧЕСКИХ СВОЙСТВ БАЛОК 


Предложено решение задачи об определении неоднородных реологических свойств балок на основе аку- 
стического метода и анализа амплитудно-частотных зависимостей. Получены операторные соотношения, 
связывающие искомые и заданные функции. Построен итерационный процесс, приведены результаты вы- 
числительных экспериментов по определению неоднородных реологических свойств материалов для раз- 
личных законов распределения искомых характеристик. 
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Введение. Материалы со сложными, в том числе неоднородными механическими свойствами 
[1-3] (полимер- и пьезокомпозиты, функционально-градиентные материалы и т.д.), широко рас- 
пространены в различных отраслях (точное машиностроение, авиастроение, измерительная 
техника, биоинженерия и т.д.). Поэтому при производстве элементов конструкций из таких ма- 
териалов и контроле качества их изготовления определение законов изменения теплофизиче- 
ских и механических свойств является весьма важной и актуальной задачей. Традиционные экс- 
периментальные методы оценки реологических свойств таких материалов в рамках гипотезы 
однородности являются достаточно грубыми, следовательно, необходима разработка альтерна- 
тивных неразрушающих методов идентификации неоднородных характеристик, позволяющая 
уточнить структуру неоднородности. 

Предложенный метод исследования основан на аппарате обратных коэффициентных за- 
дач в механике деформируемого твердого тела и позволяет восстанавливать неизвестные од- 
номерные функции по данным акустического зондирования, измеренным в некоторых точках 
исследуемого объекта [4]. Описание реологических свойств, как правило, осуществляется по 
модели стандартного вязкоупругого тела и ее обобщения на неоднородные материалы. В дан- 
ной работе рассмотрен ряд задач о восстановлении комплексной функции, характеризующей 
неоднородные вязкоупругие свойства стержня при возбуждении изгибных колебаний для раз- 
личных способов нагружения. Рассматриваемые задачи решаются с помощью итерационного 
алгоритма, основанного на аппарате интегральных уравнений Фредгольма 1-го и 2-го рода. 
Постановка задачи и построение общей схемы исследования. Исследуется задача об 
изгибных колебаниях неоднородной вязкоупругой балки длиной / при различных способах воз- 
буждения колебаний. Неоднородность механических свойств моделируется зависимостью мгно- 
венного и длительного модулей упругости от продольной координаты х. 

Уравнение установившихся колебаний в безразмерном виде после отделения временно- 
го множителя будет иметь вид [5]: 


4? _  а?иХх,к) 
а С(х, = -ки(х,ю=а, (1) 
. ике(х)+й(х) 
где С(х, 1) Ре неизвестная функция, характеризующая неоднородные реологи- 
+пк 


ческие свойства, Е(х) = Н(0)2(х), Н(х) = Н(0)й(х) — мгновенный и длительный модули 


упругости соответственно, которые могут быть как гладкими положительными функциями, 
так и иметь конечное число разрывов первого рода (что моделирует слоистые структуры). 
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ре? ЕМ 7Н(0) 
— ____ - безразмерная частота и т=и 5 
УН(о) рЕ1 





Здесь также введены параметры: к“ = 


безразмерное время релаксации. 

Для уравнения (1) рассмотрены три задачи, отличающиеся как способом нагружения, 
так и дополнительной информацией, используемой для процедуры идентификации. 

Граничные условия для различных способов нагружения можно представить в следую- 
щей форме: 

1) их0,к)=0, м(0,к)=0, Оки кю) =0, (С(Ыюи( к)’ =1 (2) 
для консольно закрепленной балки, нагруженной на свободном конце сосредоточенной силой 
при отсутствии поперечной нагрузки (4=0); 


2) и"\0,к)=0,м’(0,к)=0, С жи’, к)=Ъ (Си, к) =0 (3) 
для консольно закрепленной балки, нагруженной на свободном конце сосредоточенным момен- 
том при отсутствии поперечной нагрузки (а=0); 

3) и(0,к)=0, и"к)=0, и’(0,к)=0, и’ к)=0 (4) 
для шарнирно опертой балки, нагруженной сосредоточенной силой а = 46(х-ж). 

Соответствующие краевые задачи при заданных законах неоднородности могут быть 
решены лишь численно. Для этого был использован метод предварительного сведения к инте- 
гральному уравнению Фредгольма второго рода. Введем вспомогательную функцию у(х) сле- 

‚  а?и(х,к 
дующим образом: ъ(х,к) = бок "Ю, 
ах 

Перейдем от исходного уравнения (1) к системе уравнений второго порядка относитель- 
но функций у? и у ссоответствующими граничными условиями: 


4?и(х, к) .- у(х,к) 
4х? С(х, к) | 





и (5) 
ее =к‘и(х,к) + 4. 
Интегрируя уравнения (5) два раза по переменной х, получаем: 
Г У(5,к) 
к)= (Хх = &+Ах+А,, 
ик) | Эсе т А+ 4 
(6) 


у(х,к) = к* [© -ВихЕ, к) +4+Вх-В.. 
0 


Константы А,, 45, В,, В, находим из граничных условий (2)-(4). Подставляя выраже- 


ние для функции у(х) и меняя порядок интегрирования, получаем интегральное уравнение 
Фредгольма второго рода для определения функции и"(х,к): 


ицх,к) = к* [ (иью)К (ха Дю). (7) 
0 


В зависимости от граничных условий, ядро и правая часть имеют различный вид; для 
упругого случая этот вид представлен в работе [6]. 

Для решения интегрального уравнения (7) использован метод коллокаций [7], с помо- 
щью которого задача сведена к решению линейной алгебраической системы относительно узло- 
вых неизвестных. На основе такого подхода проведен ряд вычислительных экспериментов по 
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решению прямой задачи для различных способов нагружения (2)-(4) и различных типов неод- 
нородности; получены необходимые для решения обратной задачи амплитудно-частотные зави- 
симости. 

При постановке обратной задачи считаем известной дополнительную информацию в не- 


котором частотном диапазоне [к,,к,]: 1) о функции смещения свободного конца балки 
и(Ьк) = (к); 2) об угле поворота свободного конца и’(1, к) = Х› (к); 3) о функции смещения в 
точке приложения нагрузки х, =[0,П и(х,,к) = (к). 

Сформулированные математические задачи на основе (1)-(4) об определении комплекс- 
нозначной функции С(х1к) по описанной дополнительной информации представляют собой 
нестандартные обратные коэффициентные задачи, являются нелинейными и некорректными 
проблемами и при их решении необходимо использовать регуляризованные подходы. Нестан- 
дартность изучаемых обратных задач состоит в том, что требуется определять комплексные 
функции координаты по заданной комплексной функции частоты колебаний; существующие со- 
временные способы анализа подобных задач существенно зависят от вещественности опреде- 
ляемых функций [4]. 

Для построения комплексных операторных соотношений, которые связывают искомую и 
заданную функции, рассмотрим слабую постановку задачи. Пусть у(х,к)=и/(х,к) — гладкая 
функция, сопряженная к и„(х,к) [8]. Умножим ее на уравнение (1) и проинтегрируем получен- 
ное равенство по отрезку [0, 1]. Далее, осуществляя интегрирование по частям и используя 
граничные условия, получаем общее нелинейное операторное уравнение, связывающее функ- 
ции С(х,), их,к), Г(к), где ке[к,,к,]: 





1 
[сс п)и’(х, к)и’(х, к) к‘ихх, ких, ких = Л (к). (8) 
0 


Пусть известно некоторое начальное приближение С.(х,ж) решения уравнения (8). 
Осуществим процедуру линеаризации [9], представив функции в виде 
#(хк) = в, (х, пк) + вв (х, к),  И(хиж) = т (х, к) + й (х,), 
С(хк) = С, (х/к) + &С\ (х, к), их, к) = и (х,к)-+ 5 (х,к), 
где = — формальный параметр. 
Собирая слагаемые при одинаковых степенях = и интегрируя полученные равенства по 
частям с сохранением главного по порядку слагаемого, получаем: 





[6 (и <, киб кух = У (к)-Л(ю), кЕ[к,,к,], (9) 





где [с (хи (х, ких (Ьк)- ки (х, к)их (1, ках = Х (к). 
о 


Таким образом, сформулировано комплексное операторное соотношение (9), связываю- 
щее заданные и искомые поправки. Полученное уравнение является интегральным уравнением 
Фредгольма первого рода с гладким ядром относительно неизвестных функций первого при- 
ближения. 
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Отделяя вещественную и мнимую части для функции первого приближения комплексно- 
го модуля, получаем из интегрального уравнения (9) систему вещественных уравнений: 


[сд (®) +1 сдам бою бъкуих= Ве (© Л), 
И (10) 
[<со-поде мою озких=ний Л), 


л2к“ ск? 


1 
где введены обозначения 2,(к) = —— т, 2(кК=— у, 23(К)= 
1+ к 1+ к 


1+12к4 ° 

На основе системы операторных соотношений (10) может быть построен итерацион- 
ный процесс уточнения неизвестных функций, основанный на приближенном решении уравне- 
ний (10) с помощью метода регуляризации А. Н. Тихонова [10]. 

В соответствии с этим для различных типов граничных условий 1-й шаг итерационного 
процесса описывается следующими операторными уравнениями: 





1 я и ОИ 
[ег соль (дни! дд (®)иф бъюиь съких= Вер (®- (к), 
0 

(11) 





1 и * ЕЕ 
[ег сд-и оду м/о сыюмь боюдх= ий 9-09), 
0 


причем 2" (х) = А (5) +), = +80. 
Начальное приближение С.(х,ж) построено в классе линейных комплексных ограни- 


ченных функций, коэффициенты которых определены из условия минимума функционала не- 
вязки: 


ы. 2 
Ф= | | ,к) — Г (к) ак, (12) 


на некотором компактном множестве в четырехмерном пространстве, где и(х,,к)- функция 


смещения в заданной точке стержня при некотором законе изменения комплексного модуля. 
Процедура минимизации функционала невязки осуществлена с использованием эволюционного 
алгоритма [11]. На каждом шаге построенного процесса посредством решения интегрального 


уравнения Фредгольма второго рода было найдено новое значение и», с помощью которого 


вычислена правая часть интегрального уравнения Фредгольма первого рода и его ядро в урав- 
нениях (11). Результатом решения этого уравнения стала поправка к неизвестной функции, и с 
ее учетом был проведен следующий этап итерационного процесса; процесс останавливался при 
достижении необходимой точности в функционале (12). 

Численные результаты исследования. На основе предложенной схемы исследования для 
различных типов краевых задач по заданной амплитудно-частотной характеристике проведена 
серия вычислительных экспериментов по восстановлению функций 2(х), й(х), характеризую- 


щих законы изменения мгновенного и длительного модулей. Расчеты проведены для стержня из 
вязкоупругого материала в случае монотонных и немонотонных законов изменения реологиче- 
ских свойств [7, 11, 12]. 


Для краевой задачи (1), (2): й(х)=3-е”, ®(х) = 2 со5(п(х + 0,5)) +4, начальное при- 
ближение #,(х) = 0,6+х, 2,(х) = 1,75х + 0,75. 
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По рассчитанным амплитудно-частотным характеристикам частотный диапазон выбран 
следующим образом: [к,;к, ] = [2,5; 5,6]. 
Для достижения точности 1-10“ при минимизации функционала (12) потребовалось 42 


итерации. Значение параметра регуляризации на последнем шаге: < = 6,1625.10-° (см. рису- 
нок, а). 
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й(х) 
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1 
0 
8) 
Р(х) 
4 
2 м 
1,5 ^ - 
ке 
1 
0,5 1 
0 02 04 06 08 Хх 0 02 04 06 08 Хх 
в) 
Результат восстановления Й(х) и ©(х): 
а- для задачи (1), (2); 6 - для задачи (1), (3); в- для задачи (1), (4); —— -— исходная функция; 


- - - — начальное приближение; -+--— восстановленная функция. 
Время релаксации т =0,1 
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Для краевой задачи (1), (3): й(х) =1,4-— соз(лх-1), 2(х) = 2со8(л(х + 1,5)) +2, началь- 
ное приближение й, (х) = 0,6, 2, (х) = 3,4 . По рассчитанным амплитудно-частотным характери- 


стикам частотный диапазон выбран следующим образом: [к,;к,]=[3;5,5]. Для достижения 
точности 1.10“ при минимизации функционала (12) потребовалось 32 итерации. Значение па- 
раметра регуляризации на последнем шаге: а = 3,9812 .10-° (см. рисунок, 6). 

Для краевой задачи (1), (4): х, = 0,25, й(х) =1+ соз((х +1,5)п), 2(х) =2+2х?, началь- 

ное приближение й,(х)=1, 2,(х) =17х+2, частотный диапазон [к,;к.] = [4,4; 7,6]. Для дос- 
тижения точности 1.10“ в функционале (9) потребовалось 20 итераций. Значение параметра 
регуляризации на последнем шаге: о = 1,4946 .10-” (см. рисунок, в). 
Выводы. На основе сформулированной строгой постановки обратной задачи для дифференци- 
ального оператора четвертого порядка с помощью анализа колебаний балки решена задача об 
определении функций, характеризующих неоднородные реологические свойства материала. 
Представлены операторные соотношения, связывающие искомые и заданные функции при ана- 
лизе установившихся колебаний. Получено решение поставленной задачи на основе итерацион- 
ного процесса, сочетающего на каждом шаге решение прямой задачи и определение поправок 
на основе решения интегрального уравнения Фредгольма первого рода с гладким ядром. Пред- 
ставлены примеры определения реологических свойств для различных законов изменения ком- 
плексного модуля по длине стержня (экспоненциальных и немонотонных), для различных типов 
граничных условий. Для монотонных законов максимальная погрешность реконструкции не пре- 
вышает 5%, а для немонотонных законов максимальная погрешность реконструкции достигает 
15%, причем максимальное отклонение восстановленной функции от точного решения наблю- 
далось на незакрепленном конце в задачах (1), (2) и (1), (3) , и в точке приложения нагрузки в 
задачах (1), (4). Представленные численные результаты определения неизвестных функций по 
данным частотного зондирования, основанные на формулировке итерационного процесса, по- 
зволяют говорить о работоспособности предложенного метода в достаточно большой серии вы- 
числительных экспериментов. 

Работа выполнена при поддержке Российского фонда фундаментальных исследований 
(грант №10-01-00194-а), ФЦП «Научные и научно-педагогические кадры инновационной 
России» на 2009 - 2013 годы (госконтракт П596) и Южного математического института (Влади- 
кавказ). 
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